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AN ANALYSIS OF ENERGETIC 
SPACE RADIATION AND DOSE RATES 

SUMMARY 


The results presented in this report should be of practical value for 
preliminary engineering analysis of space hazards and planning of possible space 
missions for either near-earth orbits or deep-space probes. The dose rate 
curves shown in this report are based on the latest environmental analysis by 
NASA and the United States Air Force. The environmental models and data are 
continually updated as experimental data improves and as more sophisticated 
analyses are performed. As such modifications become available, the authors 
of this report intend to update the radiation dose rates in future reports. 


INTRODUCTION 


This study is intended to present a compendium of the environmental data 
on trapped and untrapped radiation. It considers the calculated tissue dose rates 
received from these radiations and the induced secondary radiation behind var- 
ious shield thicknesses. Energetic radiation, as used here, refers to particles 
with energies greater than about 0. 5 MeV. The results given in this report 
should be of practical value for preliminary analysis of the radiation hazard to 
man during space flight. 

The radiation environment consists of protons and electrons in the Van 
Allen belts, energetic solar radiation, and galactic cosmic radiation. The 
solar wind is not included in this study since no appreciable dose is received 
behind any nominal shielding. 

The environmental data used in this report was taken primarily from the 
works of J. I. Vette of NASA/Goddard Space Flight Center and W. R. Webber 
of the Boeing Company. 



RADIATION DOSE CALCULATIONS AND DEFINITIONS 


In this report the radiation dose rate calculations are based on the physi- 
cal property of energy deposition at a tissue point detector located at the center 
of a spherical shell of aluminum or polyethylene. The units of dose are in rads- 
tissue. By definition the rad is the physical equivalent of an energy deposition 
of 100 ergs/g ( 0. 01 J/kg) . The assumption of a tissue point detector is an over- 
simplification and always leads to a high estimate of the true dose. When skin 
dose is used, the above conditions of a point detector at the center of a spherical 
shell with no self-shielding by the target should be specified. Thus, if an astro- 
naut is placed at the point detector location he would receive considerable self- 
shielding from the rest of his body. In general, the skin doses should be about 
a factor of two lower on a man target than shown in the curves of this report. 
(This will also apply to depth doses at 5 cm or less inside the man. ) In the 
present report depth doses were derived when a sphere of tissue was placed 
about the point detector; the depth indicates the radius of this sphere. 

In this work the presented proton dose rate calculations include both 
primary protons and associated secondary particles. The secondary particle 
dose component for thin shields can be approximated fairly well by ignoring the 
attenuation of the primary flux because of inelastic collisions [ 1] . 

The units of shield thickness throughout this report are given in grams 
per square centimeter, which may be converted to centimeters by dividing this 
thickness by the material density (g/cm 3 ) . For example, aluminum p = 2.1 
g/cm 3 ; thus, 2.7 g/cm 2 of aluminum is one centimeter thick. The use of grams 
per square centimeter thickness has historical precedence in charged particle 
transport calculations and indeed is quite convenient for two reasons. First, in 
terms of this unit the penetration characteristics of charged particles are almost 
the same in all materials (hydrogen is the only exception) ; second, the mass or 
weight of a shield can be found for a simple geometry by multiplying the shield 
thickness by the surface area of the shield. 

The most important question that arises in the radiation area is the 
number of rads that a man, photographic film, or instrument can receive before 
damage is caused or detected. Ctf course, this is not always an easy question to 
answer. For example, very fast film may be damaged beyond feasible use after 
a one-rad exposure, whereas other films may tolerate 50 or more rads. Thus, 
each film type has to be investigated regarding its use and possible radiation 
exposure. Most electronic devices can tolerate very large doses before serious 
degradation of performance results. Such doses will not be encountered from 
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natural space radiations for periods of less than one year if the device is behind 
a shield of one centimeter of aluminum. Solar panels and thermal coatings on 
spacecraft may be damaged by the large surface doses they receive from low 
energy electrons and protons. 

If man is involved in space flight , his safety from radiation damage should 
be the first consideration. The following outline of human responses to an expo- 
sure in rads is given as a guide [2] . Such a memorandum , of course, cannot 
be considered as an official NASA document, but it is felt by the authors to reflect 
the extent of the radiation problem better than a table of acceptable dose levels: 

A. Gastrointestinal dose - The dose required to produce nausea 
in 10 percent of the crew = 25 rads. 

B. Skin dose - The dose required to produce erythema in 10 
percent of the crew (may produce temporary loss of hair) = 

200 rads. 

C. Hematological dose - The does required to produce a 20 
percent depression of leukocytes and/or thrombocytes = 50 
rads. 

D. Eye dose - The dose required to produce a 5 percent incidence 
of clinically significant cataracts = 80 rads. 

The writers of this memorandum suggested that the ’’Planning Operational Dose" 
should be exactly one-half the above outlined doses at the various locations in 
the body. The "Maximum Operational Dose" is that given above. The radiation 
exposure levels were selected according to the following operational criteria: 

Planning Operational Dose ( POD) : The dose which should not be 
exceeded without requiring a mission modification of some degree. The 
degree of modification will be a function of the magnitude of the excess 
dose and will be formulated by mission rules. This dose will be used 
for mission planning purposes to determine if proposed trajectories and 
time lines are acceptable. 

Maximum Operational Dose (MOD) : The dose which should not be 
exceeded without specific modification of the mission to prevent further 
radiation exposure. Such an exposure would be considered to result in a 
potentially harmful inflight response in terms of crew safety and post- 
flight response in terms of delayed radiation injury. 
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Of course, many factors enter into the determination of the accumulated 
radiation dose that a man can safely receive. For example, if he received a 
25-rad gastrointestinal dose over a two-year period the effects will be much less 
pronounced than the effect of receiving the same dose over a two-week period. 
Therefore, the above listing of human responses refer to fairly short term expo- 
sures where the human organism has not been capable of significant biological 
repairs. 

The following doses [3] are established on the assumption that the crew 
will be exposed to small increments of dose on each orbit. No allowance is made 
for pulses of radiation received at higher intensities. Table I should not be used 
for missions of longer duration than 60 days. 


TABLE I. RADIATION DOSE LIMITS FOR 30 - 60-DAY MISSIONS 


Tissue 

Depth 

POD 

MOD 

Skin 

0. 1 mm 

2. 5 rads/ day 

5 rads/day 

Eye 

3. 0 mm 

1. 25 rads/day 

2. 5 rads/day 

Bone 

Marrow 

5. 0 cm 

0 . 6 rad/ day 

1. 0 rad /day 


If one wishes to investigate missions of long duration (one or two years) 
he may assume that the body does indeed repair some of the damage; however, 
it would be presumptuous to extend the acceptable dose levels without more 
knowledge. It is conceivable that a total allowable accumulated dose may in fact 
be doubled for a mission of one or two years. Such an assumption, however, 
must embody the concept of a fairly constant or uniform radiation exposure over 
the total period. This is probably not a valid assumption for deep space flight 
since one could conceivably receive 90 percent of his allowable dose during one 
large solar proton event lasting (at most) three days. 

It is quite valid to inquire if there has been experimental verification of 
the radiation dose computational methods. The basic parameters for proton 
and electron transport in matter ( range and stopping power data) have been fre- 
quently measured experimentally and are quite valid. These basic data, however, 
are only the building blocks of the much more involved dose rate calculations. 
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The experimental verification of proton dose rate calculations, for example, 
is not easily obtained for several reasons. Laboratory proton sources are usually 
monoenergetic , whereas flare sources have continuous energy spectra. Also, 
when a theoretical computation is made for a monoenergetic source , the energy 
band is depicted as a monoenergetic line; whereas in the experiment the proton 
beam may have a 5 percent energy spread. This energy spread in the source 
may have a much greater effect on the attempt to duplicate a theoretical result 
than is commonly believed. Figure 1 illustrates the differences in proton energy 
deposition in tissue and aluminum for three different energy spreads about a 
nominal energy of 100 MeV. The source was taken as 10 8 protons/cm 2 uniformly 
distributed in the energy band 100 ± AE. The spiked curve was obtained by using an 
energy spread of only ±0. 005 MeV. This type of curve is characteristic of a 
theoretical calculation, whereas the other curves of Figure 1 are typical of 
experimental results. The only effect considered in these calculations was an 
energy spread in the source. 



FIGURE 1. MONOENERGETIC PROTON DOSE AS A FUNCTION OF DEPTH 
PENETRATED FOR VARIOUS SOURCE ENERGY SPREADS 
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In conclusion, the writers would like to warn against the indiscriminant 
use of the parametric curves presented in this study. The dose rates in general 
should be halved in order to obtain the proper dose rate on an astronaut or on 
other fairly large targets. Also, since actual spacecraft configurations are not 
spherical shells, an additional reduction in dose (as much as a factor of two) 
will be found if a realistic geometric representation is used in which shielding 
from onboard components , other astronauts , and outside structures are consid- 
ered. For very thin walled structures such as in the S-IVB workshop, such 
improvements may be very significant. For a more detailed discussion on the 
radiation transport methods, see References 1 and 4. 


THE GEOMAGNETIC FIELD AND CHARGED 
PARTICLE INTERACTIONS 


The geomagnetic field around the earth forms a "magnetic pocket" known 
as the magnetosphere. The boundary of the magnetosphere is determined by the 
solar wind which is assumed to be a radial expansion of the sun's corona. At a 
certain radial distance above the earth, the geomagnetic field energy density 
equals the energy of the solar wind and there is a breakdown of the magnetic lines 
of force. This turbulent region is the magnetospheric boundary or transition 
zone. In the magnetosphere the magnetic field dominates, while outside the 
magnetosphere, the solar wind is the controlling energy mechanism. 

The inner boundary of the transition region, called the magnetopause, 
occurs at about ten earth radii on the sunlit side of the earth, while the outer 
boundary of the transition region is in the form of a shock wave at about four- 
teen earth radii. The magnetopause around the earth is in the form of an elon- 
gated teardrop with a long tail pointed away from the sun. Figure 2 , taken from 
Ness [5] , shows a recent version of the magnetosphere. 

Perhaps the best way to represent the distribution of magnetically trapped 
particles about the earth is by using the B-L coordinate system developed by 
Carl Mcllwain [6] . The B coordinate denotes the magnetic field strength at some 
specified point in space; L is the magnetic shell parameter that labels the shell 
upon which the guiding center of the trapped particle is adiabatically confined as 
it drifts around the earth. The L coordinate is approximately constant along a 
geomagnetic field line. In a dipole field, L is constant along the field line and 
has the geometric property of being the equatorial distance from the dipole cen- 
ter to the magnetic field line. The geometry of the B-L coordinate system is 
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RESULTS OF IMP-1 MAGNETIC FIELD EXPERIMENT (11/27/63 TO (5/31/64) 



FIGURE 2. RECENT SATELLITE VERSION OF THE MAGNETOSPHERE 
BASED ON RESULTS OF IMP- 1 MAGNETIC FIELD EXPERIMENT 
(NOVEMBER 27, 1963, TO MAY 31, 1964) 

depicted in Figure 3. For a dipole field the B-L coordinates are related to the 
latitude and altitude of a point above the earth by the following: 


B 


M 

R 3 



3R \ 2 
L / 


L 


R 

cos 2 A. 


( 1 ) 


where M is the magnetic dipole moment (0. 311653 gauss R^ 3 ), R is the distance 

in earth radii units from the center of the earth, A is the geomagnetic latitude, 

and R is the radius of the earth. These relations should be used with caution 
e 

since the earth’s field cannot be depicted as a simple dipole. If the magnetic 
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FIGURE 3. THE GEOMETRY OF THE B-L COORDINATE SYSTEM 





dipole approximation is used, it should be represented by the more accurate 
tilted and translated dipole instead of the centered dipole. 


The magnetic field of the earth alters the penetration of charged particles 
to the vicinity of the earth. The theory of the allowed cone of incident charged 
particles, as developed by Stormer [7] , can be used as a basis for calculating 
the modification of an incoming energy spectrum of cosmic rays or solar flare 
protons. The theory is based on a dipole approximation of the geomagnetic field. 
The allowed cone is defined by: 

60 x 10 3 cos 4 X / 

^ ~ R 2 ( 1 + 'n/ 1 - cos^cosy j^~ 


where R and \ are defined as in equation ( 1) ; p is the magnetic rigidity (momen- 
tum/charge) in megavolts ( MV units) : and y is the half angle of the allowed 
cone about the normal to the meridian plane [8] . The angle of the allowed 
cone varies with the energy of the particle and with the position of the observa- 
tion point [9] . The modified energy spectrum inside the magnetosphere is 
found by multiplying the external spectrum by the solid angle fraction of the 
allowed cone as a function of the particle energy. 


A computer code has been developed by J. J. Wright to calculate the time 
weighted modified energy spectrum on any feasible earth orbit by varying the 
energy (angle) over the allowed cone and weighting each energy (angle) by the 
time spent at each position in orbit * . The modified energy spectra for galactic 
cosmic rays and for the November 12, 1960, solar flare will be shown in the 
following two sections. Shadowing effects of the earth were not considered in 
these cases. The magnetic rigidity in MV is related to the kinetic energy E in 
MeV of the particle as follows: 


P(E) 


\/E 2 + 2ME 

q 


(3) 


where q is the total charge, M is the particle rest mass energy (MeV) , and E 
is expressed in MeV. 

Equation 3 is plotted in Figure 4 for proton and alpha particles. For 
E < 100 MeV, p « ^2 ME/ q = 43. 5 \lE , and for E > 10 5 MeV, p « E/q. 


"I* 

For a discussion 


of methods see reference [10]. 
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KINETIC ENERGY (MeV) 


FIGURE 4. RIGIDITY AS A FUNCTION OF KINETIC ENERGY 
FOR PROTONS AND ALPHA PARTICLES 


Combining equations ( 1) , ( 2) , and ( 3) for y = 90° , and solving for the 
kinetic energy E (MeV) , the vertical cutoff energy for protons is found to be 

E = -938 + 7(938) 2 + 2,25 x ^° 8c ! os8 ^ (, 

c R* 

The so-called vertical cutoff energy signifies that particles coming from the 

zenith and having energies greater than E will intercept a given point (R,A) . 

c 

Figure 5 shows the variation of the vertical cutoff during solar active periods 
and quiet time [9] . A plot of the proton vertical cutoff energy is also shown in 
Figure 5 for various earth radii and latitudes for quiet times. The trapped 
radiation belts will be discussed further in later sections of this report. 
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PROTON CUTOFF ENERGY (M«V) 



FIGURE 5. VARIATION OF PROTON VERTICAL CUTOFF ENERGY 
DURING A SOLAR ACTIVE PERIOD AND AS A FUNCTION OF 
EARTH RADII AND LATITUDE DURING QUIET TIME 


GALACTIC COSMIC RADIATION 

The galactic cosmic radiation is composed of about 85 percent protons, 

14 percent alpha particles , and about 1 percent heavier nuclei with energies 
ranging from 10 7 to 10 19 eV, with an average energy of about 1 GeV. The free 
space proton flux at solar maximum is about 2. 5 protons/ cm 2 -sec and about 
twice this at solar minimum. 

The differential energy spectrum [11] for protons and alpha particles 
is shown in Figure 6. The galactic cosmic ray dose rate (including secondaries) 
as a function of shield thickness is also shown in Figure 6. A spacecraft having 
an average thickness of 6 g/cm 2 of aluminum would receive a total dose of 

11 


I 





SHIELD THICKNESS (g/cm z -AI.) 


b 

FIGURE 6. GALACTIC COSMIC RAY DIFFERENTIAL ENERGY SPECTRA 
AND DOSE RATES BEHIND ALUMINUM SHIELDS 


7. 2 rads/year during solar minimum. At solar maximum the dose rate could be 
reduced by a factor of two. However, if one uses the modified energy spectrum 
as computed by the technique described on page 9, the total dose from cosmic 
radiation in a 30-degree, 240-nautical mile circular orbit would be on the order 
of 0. 20 rad/year. Figure 7 shows the modified cosmic ray dose rate for protons 
plus secondaries as a function of orbital inclination on a 240-nautical mile circu- 
lar orbit for an aluminum shield thickness of 10 g/cm 2 . This graph does not 
include earth shadow effects . 


Results by some writers have 
implied that the dose rate is higher by 
a factor of two; however, insufficient 
data on their proton energy spectrum 
makes a comparison impossible. 


SOLAR COSMIC RADIATION 


The radiation intensity from 
solar flare activity follows an approx- 
imate li-year cycle with enormous 
flux variations ranging from about i0 5 
protons/cm 2 at solar minimum to 
approximately 10 9 protons/ cm 2 per 
flare at solar maximum for protons with 
energy above 30 MeV. The solar radi- 
ation intensity will also vary with dis- 
tance from the sun. The following 
solar proton data are assumed to be 
for a radial distance of one astronom- 
ical unit. 

Solar cosmic radiation has two 
main components - protons and alpha 
particles in a proton-to-alpha ratio 
varying from i to 100. There is also a minor component of heavier nuclei which 
makes up about 0. 1 percent of the total. 



FIGURE 7. TOTAL COSMIC RAY 
PROTON DOSE RATE AS A FUNCTION 
OF ORBITAL INCLINATION 



The typical intensity-time profile for various rigidities in a solar cosmic- 
ray event is shown in Figure 8. (Webber [ 12] provides detailed discussion.) 



FIGURE 8. TYPICAL INTENSITY-TIME PROFILES DURING 
A SOLAR COSMIC RAY EVENT 


The intensity follows an exponential increase to I max at t = 0 and an exponential 

decay beyond maximum intensity. The integrated intensity above a given energy 
is calculated using the characteristic rise and decay times given by 
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J(>E) = 


where I is the maximum intensity at time t = 0; t_, = rise time: and t = 
max J ’ R D 

decay time. 

The time integrated spectrum describing the flux for a given event is 
given by 


f I (>E) exp (- — \dt + f I ( >E) exp \ dt 

-i max l V o max v ‘ D ) 


(t + tj I (>E) 
R D max 


( 5 ) 


J(>p) = N 0 exp^-^-^ (6) 

where p 0 (MV) is the characteristic rigidity, and N 0 is a constant determined 
from p 0 and J(>p) . Caution must be used in applying the same spectrum down 
to low rigidity values ( <30 MeV) . 


To calculate the differential spectrum in MeV units using equation (6) 
[1] , it is sufficient to use the relationship 


= i 


+ 2ME 


(7) 


and the Jacobian 


d£ 

dE 


to obtain, 


0(E) 



p \ 

dp 

particles 

Po ) 

dE 

cm 2 - MeV - Flare 


or 


0(E) 


N n I E + M ^ ( -^E 2 + 2ME^ particles 

P ° q ' Je 2 + 2ME / 6XP V Poq / cm 2 - MeV - Flare 


where p is the rigidity (momentum/ charge) in MV units, E is kinetic energy 
(MeV) , q is the charge, and M is the rest mass energy (MeV) . Thus, for pro- 
tons, q = 1, and M = 938 MeV; for alpha rays, q = 2, and M = 3727 MeV. 


A list of the solar proton events [ 12] used in this study is shown in Table 
II. The solar proton events in Figure 9 represent the cases where the skin dose 
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TABLE II. INTEGRAL PROTON FLUX ( PROTONS/cm 2 _FLARE) AT 30 
AND 100 MeV WITH CORRESPONDING CHARACTERISTIC RIGIDITY P 0 AND N 0 


Date 

J( >30MeV) 

J(>100MeV) 

j P 0 (MV) 

N 0 

2/23/56 

1.0 x 10 9 

3. 5 x 10 8 

195 

J 

3.41 x IO 9 

8/3/56 

2.5 x 10 7 

6.0 x 10 6 

144 

1.32 X IO 8 

1/20/57 

2. 0 x 1 0 8 

7. 0 x 10 e 

61 

1. 02 x IO 10 

8/29/57 

1.2 x 10 8 

3. 0 x iO 6 

56 

8.49 x IO 9 

10/20/57 

5. 0 x 10 7 

1.0 x 10 7 

127 

3. 30 x IO 8 

3/23/58 

2.5 X 10 8 

1. 0 x 10 7 

64 

| 1. 04 x IO 10 

7/7/58 

2. 5 x 10 8 

9. 0 x 10 6 

62 

1. 18 x iO 10 

8/16/58 

4. 0 x 10 7 

1.6 x 10 e 

64 

1.67 x IO 9 

8/22/58 

7.0 X 10 7 

1.8 x 10 6 

56 

5.02 x IO 9 

8/26/58 

1. 1 x 10 8 

2. 0 x 10 e 

51 

1.21 x IO 10 

9/22/58 

6. 0 x 10® 

1.0 x 10 s 

50 

7.21 x IO 8 

5/10/59 

9.6 x 10 8 

8.5 x 10 7 

84 

1.67 x IO 10 

7/10/59 

1. 0 x 10 9 

1.4 x 10 8 

104 

1.00 x iO 10 

7/14/59 

1. 3 x 10 9 

1.0 x 10 8 

80 

2. 59 x IO 10 

7/16/59 

9. 1 x 10 8 

1. 3 x 10 8 

105 

8. 92 x IO 9 

4/1/60 

5. 0 x 10 6 

8.5 x 10 5 

116 

3. 93 X IO 7 

4/28/60 

5.0 x 10 s 

7. 0 x 10 s 

104 

5.01 x IO 7 

5/4/60 

6. 0 x 10 e 

1.2 x 10 e 

127 

3. 96 x IO 7 

5/13/60 

4. 0 x 10 6 

4. 5 x 10 5 

94 

5.09 x IO 7 

9/3/60 

3.5 x 10 7 

7.0 x 10 6 

127 

2.31x10 s 

9/26/60 

2. 0 X 10 6 

1. 2 x 10 5 

73 

5. 29 x IO 7 

11/12/60 

1.3 x 10 9 

2. 5 x 10 8 

124 

8. 98 x IO 9 

11/15/60 

7.2 x 10 8 

1.2 x 10 8 

114 

5.89 x IO 9 

11/20/60 

4. 5 x 10 7 

8.0 x 10 e 

118 

3. 44 x IO 8 

7/11/61 

3. 0 x 10 s 

2.4 x 10 5 

81 

5.77 x IO 7 

7/12/61 

4. 0 x 10 7 

1.0 x 10 e 

56 

2.83 X IO 9 

7/18/61 

3. 0 x 10 8 

4. 0 x 10 7 

102 

3. 13 x IO 9 

7/20/61 

5. 0 x 10 6 

9.0 x 10 5 

120 

3.66 x IO 7 

9/28/61 

6. 0 x 10 e 

1. 1 x 10® 

121 

4.33 X IO 7 

10/23/62 

1.2 x 10 5 

1.0 x 10 4 

83 

2. 13 x IO 7 
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a b 

FIGURE 9. TOTAL PROTON SKIN DOSE BEHIND ALUMINUM SHIELDS 
FOR SEVERAL SOLAR FLARES 


received behind 5 g/cm 2 aluminum shield is greater than about 1 rad. Figure 10 
depicts solar cosmic ray doses at the blood-forming organs. Table HI presents 
a summary from the above graphs. The proton skin dose behind polyethylene 
shields for the seven largest solar cosmic ray events is shown in Figure 11 and 
summarized in Table IV. 

Figures 9 to 11 represent free space doses due to solar flare particles. 
Because of the geomagnetic field cutoff, the doses can be reduced by one or more 
orders of magnitude. Figure 12 depicts the modified energy spectra and corres- 
ponding dose rates for the November 12, 1960, flare on a 240-nautical mile cir- 
cular orbit for various orbital inclinations. For an inclination of 45 degrees the 
dose is a constant of 0. 0003 rad/flare. These dose values do not include earth 
shadow effects. 
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TOTAL BLOOD DOSE (RADS/FLARE) 



COMPOSITE SHIELO THICKNESS (g/cm* AI./TISSUE) 



a b 

FIGURE 10. TOTAL PROTON BLOOD DOSE BEHIND ALUMINUM SHIELDS 
FOR VARIOUS SOLAR COSMIC RAY EVENTS 



TABLE III. TOTAL ESTIMATED SOLAR FLARE DOSES BY EVENT 
FOR TEN SHIELDING CONFIGURATIONS 


Shielding Configuration 

Date 

1/0* 

2/0 

5/0 

10/0 

20/0 

’ 1/5 

2/5 

5/5 

10/5 

20/5 

2/23/56 

280.00 

181. 00 

91.80 

50.20 

24.80 

64.78 

58. 00 

43.75 

30.40 

17.90 

8/3/56 

8.50 

5. 00 

2.20 

1.00 

0.40 

1.39 

1.21 

0.85 

0.53 

0.27 

1/20/57 

122.00 

43.50 

8.30 

1.80 

0.30 

3.42 

2.57 

1.23 

0.46 

0. 11 

8/29/57 


mam 

4.20 

0.80 

0. 10 

1.63 

1.20 

0.54 

0. 19 

0.04 

10/20/57 



4.10 

1.80 

0.70 

2.53 

2. 17 

1.46 

0.88 

0.41 

3/23/58 

148.00 


mm 

2.50 

0.40 

4.67 

3.55 

1.75 

0.69 

0.17 

7/7/58 

150.00 


n 

2.30 

0.40 

4.38 

3.30 

1.60 

0.61 

0. 15 

8/16/58 



1.80 

0.40 

0.10 

0.75 

0.57 

0.28 

0. 11 

0.03 

8/22/58 



2.50 

0.50 

0. 10 

0.96 

0.71 

0.32 

0. 11 

0.02 

8/26/58 



3.40 

0.50 

0. 10 

1. 19 

0.85 

0.36 

0.11 

0.02 

5/10/59 

470.00 

211. 10 

59.30 

18.30 

4.40 

30. 18 

24.28 

13.60 

6.70 

2. 10 

7/10/59 

420.00 

214. 00 

73.20 

27.40 

8.40 

41.56 

34.65 

21.76 

11.84 

4.80 

7/14/59 

650.00 

284. 50 

75.90 

22.30 

5.00 

37.56 

30. 00 

16.75 

7.80 

2.50 

7/16/59 

382.00 

194. 80 

67.20 

25.30 

7.80 

38.30 

31.98 

20. 16 

11.03 

4.50 

9/3/60 

13.00 

7.20 

2.90 

1.20 . 

0.50 

1.77 

1.52 

0. 10 

0.06 

0.03 

11/12/60 

484. 00 

269.60 

105.50 

44.90 

16.20 

64.53 

55. 12 

36.87 

21.83 

10.05 

11/15/60 

288.00 

151.90 

55.90 

22.40 

7.50 

30.04 

27.91 

18.14 

10.33 

4.49 

11/20/60 



3.60 

1.50 

0.05 

2.14 

1.82 

1.20 

0.69 

0.31 

7/12/61 



1.40 

0.30 

0.03 

0.54 

0.40 

0. 18 

0.06 

0.01 

7/18/61 



21.60 

8.00 

2.40 

12.16 

10.11 

6.30 

3.39 

1.35 


* Shielding configurations are given as X/Y where X = shielding thickness in g/cm 2 of aluminum and 
Y = shielding thickness in g/cm 2 of tissue. 










TOTAL DOSE (KAOS/FLARE) 



a b 

FIGURE 11. TOTAL PROTON SKIN DOSE BEHIND POLYETHYLENE SHIELDS 
FOR SEVEN LARGE SOLAR COSMIC RAY EVENTS 


TABLE IV. TOTAL ESTIMATED SOLAR FLARE DOSE 


BEHIND POLYETHYLENE SHIELDS 


Shield Thickness 

Date 

1* 

2 

5 

10 

15 

20 

2/23/56 

226. 8 

141. 5 

68.4 

35.7 

23. 1 

16. 5 

5/10/59 

317. 8 

134. 1 

33.6 

9. 3 

3. 8 

1. 9 

7/10/59 

303. 0 

145.7 

45.5 

15. 5 

7. 5 

4. 2 

7/14/59 

435. 0 

177. 5 

42. 1 

11. 0 

4.4 

2. 1 

7/16/59 

275. 2 

133. 0 

41. 9 

14.4 

7. 0 

4. 0 

11/12/60 

366. 1 

192.4 

69.8 

27. 5 

14.7 

9. 0 

11/15/60 

210.0 

106.0 

35.9 

13.3 

6.8 

4.0 


* Shield thickness in g/cm 2 of polyethylene . 


The November 12, 1960, proton spectral shape was somewhat controver- 
sial. The energy spectrum obtained by A. J. Masley [13] differs considerably 
from that of W. R. Webber (Table II) for energies below 100 MeV. Figure 13 
is a comparison of skin doses using both Webber's and Masley 's energy spectra. 

The differential energy spectrum used by Masley is given by 

J(E) = 1.77 x 10 13 E" 3 cm- 2 MeV -1 (30 ^ E ^ 80 MeV) 

J(E) = 9.62 x 10 16 E -5 cm -2 MeV" 1 (80 ^ E ^ 440 MeV) (8) 

J(E) = 6.63 x 10 18 E~ 5 - 4 cm -2 MeV -1 (440 =s E ^ 6600 MeV) . 

When comparing the results in the above figures, note that, using an 
aluminum shield of 6 g/cm 2 , the dose at the blood-forming organs (5 cm of tissue) 
would be approximately the same for both spectra. 

Figure 14 shows a parametric study of the dose received behind an alum- 
inum shield for various values of the characteristic rigidity, p 0 . The proton 
integral spectrum in these calculations is given by 

J(>p) = N 0 exp^-j^= 10 9 Pr c °^T ~ ( E > 30 MeV) (9) 
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DIFFERENTIAL FLUX ( protom/cm* -M«V-flor« ) 


to 

to 




a b 

FIGURE 12. MODIFIED DIFFERENTIAL ENERGY SPECTRUM AND 
CORRESPONDING DOSE FOR THE NOVEMBER 12, 1960, SOLAR FLARE 




FIGURE 13. COMPARISON OF 
USING MASLEY AND 



PROTON SKIN DOSE 
l SPECTRA 





TOTAL DOSE (RADS/FLARE) 


and the characteristic rigidity, p 0 , is given for values between 50 and 200 MV. 

Alpha particles seem to have about the same rigidity characteristics 
as protons. Figure 15 illustrates the characteristic rigidity plotted as a func- 
tion of proton-to-alpha ratio [ 14] . 



FIGURE 14. TOTAL PROTON SKIN 
DOSE BEHIND ALUMINUM SHIELDS 
FOR VARIOUS CHARACTERISTIC 
RIGIDITY VALUES 



FIGURE 15. CHARACTERISTIC 
RIGIDITY AS A FUNCTION OF 
PROTON-TO-ALPHA RATIO 


24 


J 



SKIN DOSE (RADS/FLARE) 
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To evaluate the radiation dose derived from alpha rays , the large flare 
of July 10, 1959, was used. If alpha particles have the same rigidity as protons, 
the proton-to-alpha ratio for the above flare is about 2. The alpha rays’ integral 
spectrum is given by 


J(>P) 


5 x 10 9 exp 


p \ alphas 
104/ cm 2 - Flare 


( 10 ) 


The above alpha-ray energy spectrum was used for the comparison of 
the proton and alpha-ray skin dose as a function of aluminum shield thickness 
in Figure 16. For an aluminum shield thickness of 4 g/cm 2 , the proton dose is 
a factor of 10 above the alpha dose. 



FIGURE 16. A COMPARISON OF 
PROTON AND ALPHA RAY SKIN DOSE 
AS A FUNCTION OF ALUMINUM 
SHIELD THICKNESS 


PREDICTION OF SOLAR 
PROTON EVENTS 


Because of the importance of 
solar protons in the manned space 
flight program, it seems justifiable to 
discuss the methods and status of flare 
predictions. Many statistical studies 
have been undertaken, but not too much 
reliance can be placed on these studies 
since the sample of data is rather small 
It should be pointed out that each flare 
differs from all others so that no exact 
relationship exists between observable 
features. A proper statistical analysis 
must involve a large number of events 
to permit specific statements about 
flare occurrence, duration, and inten- 
sity. 

Because of the rotation of the 
sun there exists an east-west asym- 
metry of solar proton events. For 
events occurring on the Eastern Hemi- 
sphere of the sun, the probability of 
solar protons reaching the earth is one- 
third that of events occurring on the 
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the Western Hemisphere [12] . If an event does occur on the Eastern Hemisphere 
the corresponding onset, rise, and decay times are three times greater than that 
of events on the western half, giving astronauts more time to prepare for the 
oncoming event. 

The presence and development of an active region with its associated 
sunspots and complex magnetic field is a basic part of the process which leads 
to a solar cosmic ray event. Thus , it is found that there are two aspects of 
primary importance for flare prediction [15] and warning capabilities. These 
are ( 1) the persistence of single active centers, and (2) the magnetic configura- 
tions of these active centers. Regarding the persistence of single active centers, 
Guss [16] has pointed out that a single fixed location in solar longitude produced 
most of the major events in cycle 19. During a period of over five years ( more 
than 73 rotations) , several active centers grew and died in this same local region 
on the sun. These major events included the events of February 23, 1956, July 
1959, November 1960, and July 1961. About 75 percent of the total integrated 
particle intensity above 10 MeV came from this one "hot" location. According 
to Webber [12] , over 90 percent of the output of the solar cosmic rays above 
10 MeV came from only eight major active centers during solar cycle 19. Four 
of these were associated with this one particular location. If a "hot" region exists 
and can be identified early in a solar cycle, the prediction of large events will 
probably be concerned with the study of this one region. 

According to Weddell [17] , there seems to be a linear correlation between 
the smoothed sunspot number and the number of cosmic ray events, and also with 
the integrated intensity of particles above 10 MeV. The number of annual par- 
ticle events [ 12] occurring at the next solar maximum should be about four or 
five, with an annual integrated intensity above 10 MeV of about 10 9 particles/cm 2 
based on results of solar cycle 19. According to Webber [18] the maximum sun- 
spot number should occur near the middle of 1968 and could have a value of 110 — 
130. 


The distribution of the integrated flux per event [17] as a function of the 
number of events for cycle 19 is given in Figure 17. The monthly probabilities 
of events as a function of event size over the 96-month period from 1956 through 
1963 for 68 flare events is shown in Figure 18. 

Tabulated data showed that during the solar cycles from 1942 through 
1963 (covering 76 observed events [ 17] ) , July was by far the most active month, 
with 18.4 percent; September was the second most active, with 17. 1 percent; 
August was the third most active with 13. 2 percent, and December, with no 
observed events, has apparently been the least active month. 
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Dr. J. E. Norman of the 
University of Georgia [19] in conjunc- 
tion with M. O. Burrell and J. W. Watts 
developed a stochastic model for the 
prediction of expected proton doses on 
space missions beyond the earth’s 
magnetosphere, hi this study, the 
model was based on the available solar 
flare data from cycle 19, the last com- 
plete cycle. Point skin doses inside 
spherical shells of aluminum of 2, 5, 

10, and 20 g/cm 2 thicknesses were 
calculated for each of the important 
solar flares observed during the 312 
weeks of the active period. An impor- 
tant flare was considered to be one 
which gave one rad of dose behind a 
2 g/cm 2 shield. Using this data, the 
authors developed a Monte Carlo model 
based on four assumptions. These 
assumptions were: 


FIGURE 17. DISTRIBUTION OF THE 
INTEGRATED FLUX PER EVENT 
AS A FUNCTION OF THE NUMBER 
OF EVENTS FOR SOLAR CYCLE 19 


1. The size and number 
of flares in any future 
active period were based 
on the distribution of the 
19th cycle flares. 



FIGURE 18. MONTHLY PROBABILITY OF A PROTON EVENT 
AS A FUNCTION OF EVENT SIZE 
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2. An event was assumed to consist of one, two, or three flares 
in the period of a week in the same ratios observed in the 19th 
cycle. There were 31 flares in the 19th cycle making up 25 weekly 
events consisting of 20 single flares, 4 double flares, and 1 triple 
flare. 

3. Flares making up an event were of the same size. 

4. The events were distributed randomly in time over the 6 -year 
active period. 

The individual weeks of a mission were generated by random sampling based on 
these assumptions. Twenty thousand missions were generated, and a cumula- 
tive probability distribution was formed. Figure 19 shows these distributions 
for missions of 2 and 52 weeks for the four shielding thicknesses. Figure 20 
shows dose vs shield thickness for the given cumulative probabilities. For 
example, in Figure 20 the 95 percent curve means that 95 percent of the missions 
had doses below those shown. 


TRAPPED PROTONS 


The proton environmental data used in this study were taken from the 
work of James I. Vette [20] of Goddard. The proton flux above four different 
threshold energies is shown in Figures 21 and 22. These fluxes are for circular 
orbits at different altitudes above the earth for 0-, 30-, 60- and 90-degree 
orbital inclinations respectively. The skin dose rates calculated from the AP3 
spectral data are shown in Figures 23 through 2.6 . The AP3 data were extended 
down to 40 MeV to obtain the proton dose rates for the 2 g/cm 2 cases. 

The above proton dose calculations include primary and secondary parti- 
cles, and the dose is computed for a point at the center of a spherical shell of 
indicated shield thickness. No self-shielding by astronauts is included; hence, 
doses are high by approximately a factor of two for human targets. 

Figure 27 shows an isoflux plot in the south Atlantic anomaly at an altitude 
of 400 kilometers. Figure 28 depicts the integral flux received on a 30-degree 
circular orbit at an altitude of 240 nautical miles as a function of time. 
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SKIN DOSE 0 (rod* /minion) 


b 

FIGURE 19. CUMULATIVE PROBABILITY DISTRIBUTIONS FOR 2 AND 
52 WEEK MISSIONS AS A FUNCTION OF DOSE FOR 
VARIOUS SHIELD THICKNESSES 
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a b 

FIGURE 20. DOSES RECEIVED ON 2 AND 52 WEEK MISSIONS AS A FUNCTION OF 
SHIELD THICKNESS FOR VARIOUS CUMULATIVE PROBABILITIES 


OMNIDIRECTIONAL FLUX (Protons /cm 2 -day) 



OMNIDIRECTIONAL FLUX (Protons / cm 2 -day) 



FIGURE 22. TRAPPED OMNIDIRECTIONAL PROTON FLUX ABOVE 34 AND ABOVE 50 MeV 


DOSE RATE (rods/day) 




b 


CO 
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FIGURE 23. LOW ALTITUDE TOTAL PROTON DOSE RATE AS A FUNCTION OF 
SHIELD THICKNESS AND ALTITUDE FOR 0- AND 30-DEGREE ORBITAL INCLINATIONS 




DOSE RATE (rads/doy) 




a 


b 


FIGURE 24. LOW ALTITUDE TOTAL PROTON DOSE RATE AS A FUNCTION OF 
SHIELD THICKNESS AND ALTITUDE FOR 60- AND 90-DEGREE ORBITAL INCLINATIONS 





TOTAL DOSE (Rads /day) 




TOTAL DOSE (Rads /day) 




a 
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FIGURE 26. TOTAL PROTON DOSE RATE AS A FUNCTION OF SHIELD THICKNESS 
AND ALTITUDE FOR 60- AND 90-DEGREE ORBITAL INCLINATIONS 




FIGURE 27. PROTON ISOFLUX PLOT AT AN ALTITUDE OF 400 km 
IN THE SOUTH ATLANTIC ANOMALY 



FIGURE 28. INTEGRAL PROTON FLUX RECEIVED ON A 30-DEGREE 
CIRCULAR ORBIT AT AN ALTITUDE OF 240 NAUTICAL MILES 
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TRAPPED ELECTRON AND BREMSSTRAHLUNG RADIATION 


Dr. Vette [21] has produced a projected electron environment for 
December 1968 starting with the August 1964 AE2 electron environment. He 
decreased the number of electrons in the inner belt according to the measurements 
of Bostron, et al. [22] , and increased the number of electrons in the outer belt 
to correspond to solar maximum conditions. Figure 29 shows the 1968 data as 
compared to the 1964 electron data. Figures 30 to 33 show the electron and 
corresponding bremsstrahlung dose rates for the 1968 electron environment. 

( Bremsstrahlung is the electromagnetic radiation produced by the acceleration 
and deceleration of electrons in nuclear electric fields. ) 

The 1964 electron environmental data and the related electron and brems- 
strahlung dose rates have been placed in the appendix. The 1968 data are now 
considered to be the best available environmental model for planning purposes. 

The development of a new electron environmental model is being pursued by 
Dr. Vette at Goddard, but may not be forthcoming until the effects of the present 
solar cycle are better known. 


A typical bremsstrahlung dose transmission curve is shown in Figure 34 
for isotropically incident electrons with an energy spectrum given by 


<M>E) 


10 13 exp( -2E) 


electrons 
cm 2 - day 


(ID 


The bremsstrahlung dose is not strongly dependent on aluminum thickness, as 
shown in Figure 34. However, the bremsstrahlung is primarily of low energy 
and its intensity can be greatly reduced by adding a thin layer of lead on the 
inside of the shield. For shields thicker than about 2 g/cm 2 the bremsstrahlung 
dose is the most important radiation hazard associated with electrons. The 
isotropic incidence assumption is more valid for electrons and associated brems- 
strahlung dose calculations than a simple straight ahead approximation. 

The electron and bremsstrahlung dose rates as a function of shield thick- 
ness are shown in Figure 35 for a model electron integral spectrum given by 

<M> E) = N Q expf-|j- j (12) 

° 'o' 

with 0(> 0. 5) = lO 10 electrons/cm 2 - day and the values of E are chosen between 
0. 25 and 1. 25 MeV. ° 
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OMNIDIRECTIONAL FLUX (ELECTRON/cm -DAY) OMNIDIRECTIONAL FLUX (ELECTRON/cm- DAY) 


r 



ALTITUDE (I0 3 nmi) ALTITUDE (I0 3 n mi) 

COMPARISON OF ORBITAL INTEGRATIONS COMPARISON OF ORBITAL INTEGRATIONS 
AT 0° INCLINATION. AT 30° INCLINATION. 


AE2 MAP AUGUST 1964 

PROJECTED 1968 MAP DECEMBER 1968 



0 3 6 9 12 15 18 0 3 6 9 12 15 18 


ALTITUDE (I0 3 nmi) ALTITUDE (I0 3 nmi) 

COMPARISON OF ORBITAL INTEGRATIONS COMPARISON OF ORBITAL INTEGRATIONS 
AT 60° INCLINATION. AT 90° INCLINATION. 

FIGURE 29. PROJECTED ELECTRON ENVIRONMENT FOR DECEMBER 1968 
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ELECTRON SKIN DOSE RATE (rods/doy) 




a b 

FIGURE 30. PROJECTED ELECTRON AND BREMSSTRAHLUNG DOSE RATE 
AS A FUNCTION OF ALTITUDE AND SHIELD THICKNESS FOR 0 -DEGREE ORBITAL INCLINATION 


ELECTRON SKIN DOSE RATE (rad*/day) 




a b 

FIGURE 31. PROJECTED ELECTRON AND BREMSSTRAHLUNG. DOSE RATE 
AS A FUNCTION OF ALTITUDE AND SHIELD THICKNESS FOR 30-DEGREE ORBITAL INCLINATION 


ELECTRON SKIN DOSE RATE (rodi/doy) 



FIGURE 32. PROJECTED ELECTRON AND BREM SSTRAHLUNG DOSE RATE 
AS A FUNCTION OF ALTITUDE AND SHIELD THICKNESS FOR 60-DEGREE ORBITAL INCLINATION 
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FIGURE 33. PROJECTED ELECTRON AND BREM SSTRAHLUNG DOSE RATE 
AS A FUNCTION OF ALTITUDE AND SHIELD THICKNESS FOR 90-DEGREE ORBITAL INCLINATION 



SKIN DOSE (RADS/DAY) 



0 1 2 3 4 5 6 7 

Al. SHIELD DEPTH ( g /cm 2 ) 


FIGURE 34. ELECTRON AND BREMSSTRAHLUNG DOSE RATE 
AS A FUNCTION OF SHIELD THICKNESS USING A 
TYPICAL ELECTRON INTEGRAL SPECTRUM 


44 



MODEL ELECTRON DOSE (RADS/DAY) 


10 ' 


K ) 3 


5 j 

MODEL SPECTRUM 




Ol 


FIGURE 35. ELECTRON AND BREMSSTRAHLUNG DOSE RATE AS A FUNCTION 
OF SHIELD THICKNESS FOR SEVERAL CHARACTERISTIC E q VALUES 
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SYNCHRONOUS ORBIT RADIATION 


The galactic cosmic rays, solar flare protons, and trapped electrons are 
the important dose sources at the synchronous orbit altitude. The dose rates for 
galactic cosmic rays are shown in Figure 6 of a previous section. The free space 
dose rate from solar proton events were shown in Figures 9 to 11. Since the 
vertical cutoff for protons at the synchronous altitude is about 60 MeV, the solar 
and galactic cosmic ray dose rates will be only slightly modified. 

Figure 36 depicts the synchronous orbit electron integral flux as a function 
of energy and local time as taken from Vette [21] . This electron integral flux 
is given by the equation 

_ p /Tp 

, . „ „ . . _ . . .. . „ . -b x o N(t) electrons . 

<M>E, B/B q , t) = C • A(t) (B/B q ) e E cm 5 sec" 13) 


where 

C = 9 x 10 7 at solar minimum ( 1/2 this at maximum) 

b = 0.625 

E = 0. 215 MeV 
o 

This equation is valid for energies greater than 0. 01 MeV. The values of A(t) 
and N(t) for various local times are given in Reference 21. The values of the 
above constants at 12:00 noon are A(t) = 1 .00 and N(t) = 0. 000. 

When the above equation is used, the synchronous altitude electron and 
bremsstrahlung dose rate as a function of local time and shield thickness would 
be as shown in Figure 37. The dose rates obtained by integrating over the curves 
in Figure 37 are shown as curve A of Figure 38. 

The time average electron integral flux at the equatorial synchronous 
altitude is given by Vette [21] as 

<M>E) = 5. 2 x 10 7 exp ( -5E) (14) 

cm sec 

The electron dose rate when this average electron spectrum is used is shown as 
curve B in Figure 38. The bremsstrahlung dose rate is shown as curve C in 
this figure when the above average flux is used. 
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FIGURE 37. SYNCHRONOUS ORBIT ELECTRON AND BREMSSTRAHLUNG DOSE RATE 
AS A FUNCTION OF SHIELD THICKNESS AND LOCAL TIME 





ELECTRON DOSE RATE (rads/day) 



SHIELD THICKNESS (g/cm 2 -AI.) 


FIGURE 38. SYNCHRONOUS ORBIT ELECTRON AND BREMSSTRAHLUNG 
DOSE RATE USING AN AVERAGE ELECTRON SPECTRUM 



CONCLUSIONS 


The results presented in this report should be of practical value for 
preliminary engineering analysis of space hazards and planning of possible 
space missions for either near-earth orbits or deep-space probes. The dose 
rate curves shown in this report are based on the latest available environmental 
analysis by NASA and the United States Air Force. The environmental models 
and data will be continually updated as experimental data improve and as more 
sophisticated analyses are performed. As such modifications become available, 
the authors of this report intend to update the radiation dose rates in future 
reports. This is pertinent to the synchronous orbit radiation environment that 
is now being examined in greater detail by experimental space probes and 
theoretical analysis. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Huntsville, Alabama, November 2i, 1967 
124-09-11-00-62 
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APPENDIX 


The 1964 electron environmental data, now considered to be out of date, 
have been placed in the appendix. Figure 39 shows the August 1964 electron 
flux data as a function of altitude and orbital inclination. The electron and 
bremsstrahlung dose rates that are shown in Figures 40 and 43 make use of 
these data. 



FIGURE 39. TRAPPED ELECTRON FLUX FOR CIRCULAR ORBITS 
AS A FUNCTION OF ALTITUDE AND ORBITAL INCLINATION 
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FIGURE 40. ELECTRON AND BREMSSTRAHLUNG DOSE RATE AS A FUNCTION OF 
ALTITUDE AND SHIELD THICKNESS FOR 0-DEGREE ORBITAL INCLINATION 
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FIGURE 41. ELECTRON AND BREMSSTRAHLUNG DOSE RATE AS A FUNCTION OF 
ALTITUDE AND SHIELD THICKNESS FOR 30-DEGREE ORBITAL INCLINATION 
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FIGURE 43. ELECTRON AND BREMSSTRAHLUNG DOSE RATE AS A FUNCTION OF 
ALTITUDE AND SHIELD THICKNESS FOR 90-DEGREE ORBITAL INCLINATION 
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